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Abstract 

The Hamiltonian of the ./V-particle Calogero model can be ex- 
pressed in terms of generators of a Lie algebra for a definite class 
of representations. Maintaining this Lie algebra, its representations, 
and the flatness of the Riemannian metric belonging to the second 
order differential operator, the set of all possible quadratic Lie alge- 
bra forms is investigated. For N = 3 and N = 4 such forms are 
constructed explicitly and shown to correspond to exactly solvable 
Sutherland models. The results can be carried over easily to all N. 
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1 Introduction: The Calogero model 



The Calogero model is a quantum mechanical system of N particles mov- 
ing on a line. It is completely integrable and possesses both fermionic (totally 
antisymmetric) and bosonic (totally symmetric) solutions. For a review see 
0. In [Q] it has been shown that it can be reduced to a representation theory 
problem of Lie algebras, so that both eigenvalues and eigenfunctions can be 
extracted from finite dimensional representations of these algebras. 

The idea that exact solvability and Lie algebraic representation might be 
connected was presented in ||, which dealt with exactly and quasiexactly 
solvable elementary models. Their Hamiltonians were rewritten as quadratic 
polynomials in Lie algebra differential operators. In this article it was also 
discussed first whether this procedure could be inverted: Find new exactly 
solvable models from Lie algebras of differential operators. Usually this leads 
to kinetic energy terms which are Laplace-Beltrami operators on a curved 
space. So the task arises to construct quadratic polynomials from the Lie 
algebras with the constraint that the curvature tensor vanishes identically on 
some dynamically accessible domain ("flat quadratic Lie algebraic forms"). 
In order to maintain the notion of exact integrability, however, the nontrivial 
but curvature-free Laplace-Beltrami operator should be transformed into a 
standard Laplace operator by an explicitly known diffeomorphism. This has 
to be found from nonlinear partial differential equations. In the neighborhood 
of the Calogero models for iV = 3 and N = 4 such approach is feasible as 
we shall show in this article. The solutions we give depend on a deformation 
parameter s ranging over the full real axis, the case s = belonging to the 
Calogero models themselves. 

In the last year two further articles have appeared dealing with the alge- 
braization of solvable models |5], || . We hope that these results can be used to 
simplify the inverse issue, the construction of new solvable dynamical models. 

Let us now outline the algebraization of the standard (sl(N)) Calogero 
model, since we use it in the sequel. The Hamiltonian is 

H = + 1 -l-A + ^]+j: 7 ^ JL rY 2 (1-1) 



where x G Rw and 



g is a coupling constant 



so that 



^ i<j \ X i x j, 



N 

* 2 = 5>* (1.2) 

i=l 



g = v{y -1)>~ (1.3) 



V = - ± (g+ -)a 
2 vy A ! 
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is an integer. If v is odd (even) we find totally antisymmetric (symmetric) 
solutions. 

We factorize the solutions as 

${x) = V{x) u e~^ x2 P{x) (1.4) 

where V(x) is the Vandermonde determinant 

V(x) = - Xj) (1.5) 

and P(x) is a polynomial which is totally symmetric in {xi}± . If ip is a 
solution of 

#V = £ty (1.6) 

then P has to solve 

hP = eP (1.7) 

with 

*--*+*S*5r"§^(=r4) (L8) 

and 

e = 2£ - iVa; - z/iV(iV - l)u; (1.9) 
Next we introduce elementary symmetric functions (following |j|) 

N N 

11(1 + atf) = £ (7„(x)r (a = 1) (1.10) 

1=1 71=0 

For each element g G (the symmetric group of N elements) we have a 
"sector" 



Eg ^■^'11 



(1,2,3,... AT) _> (ix.ia,...^)} 

ClR^v (1.11) 

so that the map 

o- : * - K}£U (1.12) 
is a diffeomorphism. This follows from the fact that the Jacobian 

M* = ^ (1.13) 

has 

r Ni 

detM = V(x)(-l) [ - ] (1.14) 
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where V(x) is the Vandermonde determinant. Its square is a polynomial 
v (<7i, . . . , o"at). The domain D + C IRat where v > is the image u(E g ) for 
each g E Sjy. 

Next we eliminate the centre-of-mass motion by the substitution 

1 

x i ^>y i =x i - — oi (1.15) 



Since 
implies 

we have 
for 



N N 

n(i+^) = i+E^ n ( li6 ) 

i=l n=2 



~ r iV i 

detA^ = 1/(x)(-l) [ - ] (1.18) 



^ = f ^ ^ ^ Tjv l n 19) 



Therefore we can use 

{(Ji,r 2 , . . . ,r n } 

as coordinates on each E g . 

On Eg we have the Laplacian (see 0) 



Q2 N Q2 N q 

A=w *? + £,^ + S^ (L20) 



with 



+ E (* - 3 - 20r i+ i_i7*_ { _i (1.21) 

Z>max(l,fc— j) 

B 4 = -i(iV - z + 2)(iV - i + l)r,_ 2 (1.22) 

where we insert 

r = l, Ti = for all z £ {0, 2, 3, . . . , N} (1.23) 
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For h fll.81) we have then 



qi d N d 2 



Serf " 9ai ^ 2 dTjdn 
^ d / 1 \ ^ d 

+**Ztow i + for + ") g (Ar - 1 + 2)(/v - 1 + 1)Ti - 2 ^ (L24) 

Next we introduce the Lie algebra of operators 

^=g^, J ij = Ti-^> i,j e{2,3,...,N} (1.25) 

This algebra has the structure 

gl(N — l)[s]£/v-i (t N -x abelian) (1.26) 
We represent it in the polynomial spaces V n 

U m = spanjrf ,r 3 ^,...,rr,E^ = m } (1-27) 

dimV n =( n + N ~ 1 ) (1.28) 



n 

so that 

V^o C V 1 C y 2 C . . . 

forms a flag of polynomial spaces. The algebra (|1.25|) can be obtained from a 
maximally degenerate representation of sl(N) (Young tableau with only one 
row). 

Finally we separate off the centre-of-momentum motion by 

+ 2 Lua — 



h = -N— + 2ua 1 — + h rel (1.29) 



and obtain the relative motion Hamiltonian h re \ as a quadratic polynomial 
in the eveloping algebra of (|1.25|) , ( |1.26|) (see |£|). 



2 The inverse issue 

We consider a Lie algebra spanned by a basis { Ja}f f being respresented in a 
flag of polynomial spaces 

V C Vi C V 2 C . . . 
rik = dim\4 < oo 
1 = n < rii < n 2 < . . . (2.1) 
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and such that 

J a V n C V n (2.2) 
Then the spectrum and the eigenvectors of the form 

M M 

^ C a pJ a Jf3 + 2_j G a Ja (2-3) 
a,/3=l a=l 

with constant coefficients {C^, G q } can be calculated from representation 
theory. 

The quadratic term in ( |2.3| ) defines a Laplace-Beltrami operator which 
we want to relate with a standard Laplacian from the kinetic energy term 
in a Schrodinger equation. For this purpose it must have zero curvature. So 
the first task is to find the flat quadratic Lie algebraic forms 

EVa^ (2.4) 

a,/3 

This is a highly non-trivial problem. 

We assume that a curvature free quadratic form 

^CtyjaJp (2.5) 

(one point in quadratic Lie algebraic form space) is known. For t G Mm 
define 

U(t) =exp{]Tt a J a } (2.6) 



and its action 



Then 



U(t)J a U(t)- 1 = Y J T a p(t)J (3 (2.7) 

& 



a/3 J 

= Y.C^{t)J a J p (2.8) 



with 



C aP (t)=J2C ( $T ia (t)T 5p (t) (2.9) 

7,5 

defines an M-dimensional orbit space of flat quadratic Lie algebraic forms. 
In the case ( |Q5l) , ([OBD 



^ = Z7(t)r fc E7(t)- 1 = e fc (t;r) (2.10) 
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is an affine linear mapping. The flat quadratic Lie algebraic forms we are 
searching for are determined up to such Lie algebraic automorphisms. 

(0)- 
xf3. 



Assume now that {C^l} is contained in a curve {C Q/ g(s)} of flat quadratic 



Lie algebraic forms for s = 0. Since {C^l} is assumed to be flat, there exists 
a diffeomorphism of some dynamically accessible domain 

7* = ie{l,2,3,...r} (2.11) 

so that 

r q2 

C^p J a Jj3 = tt~9 + lower order diff. operator (2.12) 

a,/3 i=l ^ 

Then we need a 1-parametric set of diffeomorphisms 

7* = fi(s;r') 

i G {1,2,3,. ..r} (2.13) 

so that 

CqpJgj^J^T) = C^l J a (r') Jff (t 1 ) + lower order diff. operator 

a,/3 a,j3 

(2.14) 

Both types of diffeomorphisms (|2.11|) , ( |2.13|) ought to be known explicitly, 
in order to maintain the concept of "exact solvability". 

In the following sections we shall give a curve of forms for the Calogero 
point {C^g} for iV = 3 and N = A satisfying these requirements. 



3 The N = 3 Calogero model and the flatness 
constraint 

In the case N = 3 we obtain the following Lie algebraic operators from (|1.25|) 
in lexicographic order 

d d _ d 

J\ — TJ — j J 2 — Tj — ) J 3 — T 2T\ — 

dr 2 ot 3 ot 2 

d d d 

Ja = t 2 t^, -h = n—, J & = r 3 — (3.1) 

OT 3 OT2 OT 3 

and from ( |L20|) , (|L"2T|) (then (0) belongs to -h Q) and +A ( |1~20D ) 

C{? + C$g = -3 (3.2) 

= + 2 - (3.3) 

CS } = "I (3-4) 
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all other coefficients are vanishing 

r (0) , ^(0) _ r (o) r (o) _ n 
°14 °23 — °36 "T °45 — u 



(3.5) 



^(0) _ ^(0) _ ^(0) _ ^(0) _ ^(0) _ ^(0) _ 
Wl — ( - y 12 — ^15 — <~ y 22 — °24 — ( - y 26 _ 
^(0) _ ^(0) _ ^(0) _ ^(0) _ ^(0) _ ^(0) _ ^(0) _ n ( n R , 

°33 _ °34 — °35 _ °46 _ °55 _ °56 _ °66 _ u \°-°) 

In order to find the set {C a/ g(s)} corresponding to a flat quadratic Lie 
algebraic form we can calculate the curvature tensor, each component of 
which is a rational function of the {C a /3{s)}. Equating a numerator to zero 
means that in turn all coefficients of the monomials t^t™ vanish. This leads 
to an enormous number of equations of high polynomial degree in the 21 
variables {C aj g(s)}. However, after we discovered a solution this way, we 
developed a simpler algorithm whose variable set is considerably reduced. 

First we take account of the irrelevant automorphisms (|2.10| ) and postu- 
late that ( [2.131 ) takes the form 

r i = r; + 0( S ) ie{2,...,N} (3.7) 

in general. Second we make an ansatz 

N d d 

a,/3 a,b=2 OTa OTb 

+lower order differential operator (3-8) 
where each <?<^ (s; r) is at most of degree two in {r n }. 

9ab( S ; T )= Y, lab^T^j (3.9) 

i,je{0,2,3,...,jV} 



(ro = 1) 



For s = we postulate 



9ab(0;r) = (g^)-b 1 (r) (3.10) 
and in particular for N = 3 from ( |3.2j )-( ]3"lj| ) 

<W = ( 111 ^ ) (3.1D 
Third we introduce the concepts of a "dimension" A n for each r n 

[T n ] = A n (3.12) 



S 



and of a "signature" 

sign(r 2ri ) = +1, sign(r 2n+ i) = -1 (3.13) 

Then the whole expression ( |3.8| ) is (from ( |3.11| )) required to be of dimension 
A^ 1 and of sign +1. The existence of independent dimensions (|3.12| ) for 
each r n is connected with the fact that (|2.6| ), ( |2.10| ) permits separate dilation 
automorphisms 

U n (t) = exp j*r n ^-J (3.14) 

The parameter s has naturally a dimension (see below) 

s:[s} = \ 2 1 (3.16) 
But besides s we need other dimensional quantities for n G {3, . . . , N} 



w n : [w n ] = -p- (3.17) 



Aj 
A,, 

The set 

{s,w n }% =3 

contains the "scaling parameters". At the end it will turn out that the 
diffeomorphism ( |2.13j ) can solely be expressed by functions of scale-invariant 
variables 

T} n = w n szT n , n e {3, . . . ,N\ (3.18) 

Using these concepts the most general ansatz for (|3.9| ) is for N = 3 

u x ( 2r 2 + sr| + a 2 w|s 2 r 3 2 , 3r 3 + a lS r 2 T 3 \ 
9ab U ' T) ~ " i 3r 3 + a lS r 2 r 3 , -|^ 2 r| + a 3 srl ) ^ 

Obviously neither s nor w 3 can be determined from the flatness constraints. 
In order to determine the three numerical parameters a\ , a 2 , a 3 we may set 

s = 1, w 3 = 1 (3.20) 

respectively we must obtain ([3.1 1|) for 

s = 0, w 3 = 1 (3.21) 

Then it results from flatness 

4 1 

Ol = o' ° 2 = fl 3 = 1 ( 3 - 22 ) 

3 o 
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4 The basic diffeomorphism 

The basic diffeomorphism ( |2.13p is obtained from the system of nonlinear 
differential equations 

J- 1 (^)=fi r (s W )"W (4- 1 ) 
where by derivation of ( |2.13 ) 

V ab (T') = fj(r') (4.2) 

a 

We use the scale invariant variables ( |3.18| ) with 

V3 = V ( 4 -3) 
Then the Laplace-Beltrami operators appearing in ( 2.14| ) are 

d d 



^' lh} ) I ^° differentiation ofr ) ( 44 ) 



with 



(4.6) 

The diffeomorphism ( [2.13|) involves two functions 



V = G{?,rf) (4.7) 



that must be determined from 

T^,r ] ) = H T T^\C,V r )H (4.8) 

where 

/ OF 8G \ 

H=( d £ 1 (4.9) 

V drj' drj' / 

Imposing the boundary conditions ( |3.7| ) the solution of ( |4.8| ) is unique and 
can be given in series form 

oo oo 2 n+3m ~ 1 (v -\- 2m — IV 

F(z,ri)= y y(-i) m 4 — } :x n v 2m (4.10) 

/; ^ ^ 1 n\(2m)\(2n + 6m - 1)P ' V ; 

n=0 
(n+m>l) 
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oo oo cyn+im+lr _i_ OmV 

g&v) = EE(-ir ,, 9 -r^ \T\o^ 2m+1 ( 4 - n ) 

^o^o n!(2m + l)!(2n + 6m + 2)! 

Both functions are entire analytic. The Jacobian H ( [4.9| ) is, however, not 
throughout invertible on the whole £,?7-plane. 

A summed up form for F and G is obtained after insertion of the functions 
(see CT), CT) ) 



namely 



£ = s(yxy 2 + 2/22/3 + VsVi) (4.12) 

3 

v = s^ymyz (4.13) 



ill 

S \ 2 - , / S \ 2 - r / S \ 2 



cycl. perm, of 1,2,3 



2 E sin{(jj W}™{UJ fc> «*{(§) ^}( 4 - 14 ) 



i 

S \ 2 . r/S\2 , . r / S \ 2 



G = 8^ sin {(-j J/Jsiniy y 2 }sm{f-) y 3 } (4.15) 



From ( |4.10| ) , ( [4.11| ) one can easily get the series expansion for det H 



detH(^ V )= £ h nm C V 2m (4.16) 

n,m=0 
(/KM = 1) 

but a closed form for h nm has not been obtained. Nevertheless we were able 
to sum this series as follows. 

The square of the Vandermonde determinant gives (at w% = 1) 

D = 4£ 3 + 27r] 2 = —s 3 V(y) 2 (4.17) 

so that s > (s < 0) implies D < (D > 0) respectively. The separatrix 

D = 

separates a domain 

.D < 0: here det H = can be fulfilled on an infinite number of zero-trajectories 
C n , ne {1,2,3...}; 

from the domain 

D > 0: where det H > 0. 

The curves C n are given by 

/ 3 £ \ 2 D 

C n- U + o-h +^T = ( 4 - 18 ) 
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Figure 1: The zerotrajectories C n and the separatrix in (£, 7/)-plane 



They were first discovered by numerical computations (Fig. 1). 

It is obviously suggesting to attempt a WeierstraB product ansatz 

/ \ 2 

dettf(£,77)=n 



n=l 



3 f 



+ 



D 



8n 6 7r 6 



That this is in fact correct can be verified in the following fashion. 
Consider the polynomial 

x 3 + ux + v = 

with its factorized form 

3 

Y[(x - ctj) = 
1=1 

so the elementary symmetric polynomials of {a}f =1 arc 
(J\{ol) = 0, 02(a) = u, 03(a) = — v 
By symmetric evaluation of the product we have 



n 

n=— 00 
(symm. 
n^O) 



n7r 



n [1 

n=l \ 



a. 



n 2 7r 2 



(4.19) 



(4.20) 
(4.21) 

(4.22) 
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Correspondingly we get 
det H(Z,7i) 



n 



n=—oo 
(symm. 
n^O) 
3 



2 n 2 7T 2 



-£> 1 

8 n 3 ir 3 



n 

i=i 



sin a,- 



ft; 



where 



it = v 



(4.24) 



(4.25) 



We can expand ( |4.24| ) into a power series of fti, ft2, ft3- Symmetry makes sure 
that this series can be reformulated in terms of 

(Ti(ft), cr 2 (ft), cr 3 (ft)- 



Using then Q4.22j ) with ( |4.25| ) we obtain a series expansion as in (|4.16|) . By 
a numerical comparison we can prove their identity. 
Finally we consider the separatrix 



D = implying (say) a 3 = 0, ft 2 = — fti 
This yields by (|422|), ( ggg ) 



cr 2 (a) 



fti 



Then 



det #(£,t?)|d=o 
implying zeros of det H at -D = and 



sin fti 
fti 



-e = (nvr) 2 , n E 7Z - {0} 



(4.26) 



(4.27) 



(4.28) 



(4.29) 



At this point the zero trajectory C n touches the separatrix. 

The representation of the function det H (£,17) simplifies if we return to 
the relative coordinate variables yi,y2,V3- The cubic equation (|4.20|) suggests 
the identification 



ft; 



V2 



s2 (Vj -Vk) 



(4.30) 
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k cyclic) 

Then we find easily 



o"i(a) = 

3 3 
cr 2 (a) = -sa 2 (y) = -f 



aM = -sim = ^=R (4.3!) 

so that ( |4.22j ), ( |4.25| ) are satisfied. The Vandermonde determinant V(y) was 
defined in ( |1.5| ). In y-sp&ce the separatrix consists of three straight lines 

Hi ~ Vj = 0, all pairs i ^ j (4.32) 

cutting the plane into six sectors. They correspond to the six elements of the 
symmetric group S3 (see ( |1.11| )-( |TTT4"| )). From ( ^4.24] ), ( |4.3U| ) we see that the 



inverse image of a zero trajectory C n consists of six straight lines as well 

~^ s ^(yi ~ %) = n7r ' n G IN ' a11 p airs (*>^) * ^ j ( 4 - 33 ) 

Three of these lines intersect the same sector and these lines intersect each 
other on the boundary of the sector, i.e. on the separatrix. In Fig. 2 we 
present the resulting image for the identity element of S3 as sector, projected 
on the (1/2 > 2/3) -plane. 



5 The spectrum 

On the linear space ( 1.27 ) 



V n =®U n (5.1) 



m=0 



we consider the quadratic operator 



in mv)(i) (5.2) 



017 



\ <9£ ' drj j 

(differentiation of V included) with V the matrix ( |4.5|) . This operator decom- 
poses in two parts A + B so that for all m 

AU m C U m , BU m+1 C U m (5.3) 

Now we span U m by monomials of degree m 

m 

a^Cv m ~ r (5-4) 

r=0 
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Figure 2: The sectorial in (3/2, 2/3) -plane with zerotrajectories 1 < n < 5 
then the map .4 is 



4 m) ^ < {m) = [m 2 + \m + |r(m - r)]a< m) 



— (m — r + 2)(m — r + l)c4™ 2 
3 



-i(r + 2)(r + l)aftJ (5.5) 

4 m) ^ af^ 1 ) = (r + l)(6m - 4r - l)a<+{ (5.6) 

From (|5.5|), (|5.6|) we recognize that all spaces U m decompose into direct sums 
of even and odd parts (the power r of £ is even or odd, respectively) 

U m = U<p © U£) (5.7) 

U^ = s pa n{e r V m - 2r M (5-8) 

C/H = spanie 2 ^ 1 ^- 2 ^ 1 }!^ J (5.9) 

so that 

•AU^ ^ C ^ 

BC/W c (5.10) 
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Correspondingly the set of eigenvalues A^ of A on U m decomposes 

A ra = AL +) UAH (5.11) 

The eigenvalues on each subspace Uffl are non-degenerate for fixed m. 
For the even part we find 

A« = {him 2 - (4r - 10)m + 4(r - l) 2 ]; 

o 



r6 {l,2,...,[-]+l}} 
For increasing r the eigenvalues A^ 71 -' decrease. Moreover we find 

A^ = A^ — {m 2 + 2m}, m even 

where m 2 + 2m is the smallest eigenvalue of Affl . 
On (|5.1|) A + B has the matrix form 



(5.12) 



(5.13) 



/ A 











B n , 




-1) 










B(n- 


1) 


A(n- 


-2) 









B {n . 


-2) 





l(n~3) 



V 



o 

£ 2 

o Bi a y 



(5.14) 



where A m {B m ) is the restriction of .A(j6) to C/ m . Assume we want to calculate 
an eigenvector in V n which is not in V m (m < n). If 



A m n A m i 



for m ^ m! 



(5.15) 



then 



Ae„ = A( n) e n , e n eU n (5.16) 
implies for the component e m of the same eigenvector in U m , m < n that 

e m = -{Am - A( m) E)- 1 S (m+1) e (m+1) (5.17) 

However, (|5.15|) is not correct. Instead we found an infinite sequence of 
counter examples. The smallest ones are 

- 1026 



A< 14) 


— A 8 


868 
3 


Af 7) 


— A 15 


Af> 


,(20) 
— A 9 


1336 
3 


A 6 30) 


— A 15 


,(22) 
A 3 


\(24) 
— A 10 


= 636 


\t ] 


,(32) 
— A 12 


a! 22) 


x(24) 
— A 6 


2068 
3 




= Af } 


Ai 26) 


,(28) 
— A ll 


_ 2584 
3 




,(35) 
— A 16 



3280 
3 

3364 
3 

3664 
3 

3910 
3 



(5.18) 
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Two eigenvalues can only coincide if mi — m 2 is even and m^rns are large. 
Consequently for such coinciding pair with mi < m^ri < r 2 the inverse 



(An: - X^E)- 1 

does not exist. But in addition 

^mi+l6mi+l 

^ (5.19) 



since B m ( |5.6|) is diagonal and its elements are all non-zero. 

Thus the eigenvectors in V m2 belonging to A^ 2 - 1 are missing. Any smooth 
function in £ and r\ which can be approximated uniformly over compact do- 
mains by polynomials due to the Weierstrafi approximation theorem cannot 
be approximated by the eigenpolynomials of the operator ( |5.2|) . 

This incompleteness can be cured by adding a linear term in the Lie 
algebra to (p.2|), e.g. 

+ 7 (^f + %)' 7 ^ (5 - 20) 

In fact we shall see in the subsequent section that such term is essential 
for self-adjointness of the Schrodinger operator. 

6 The Schrodinger operator 

We start from ( |3.19| ) and consider the operator 

a,b a ° a u 1 a 

where 7 € 1R is arbitrary and the functions r a (r) are linear in {r a } and left 
otherwise open for a while. The spectrum of Q can be calculated with the 
method of Section 5. Then we may have eigenfunctions if and eigenvalues e 

Q<p = e<p (6.2) 

By an appropriate choice of a gauge function \ 

(p = e x il) (6.3) 

we want to transform Q into a standard Schrodinger operator in r-space 

- + Wip = eip (6.4) 

Here A is the Laplace-Beltrami operator 

A =^5£^£ (6 - 5) 
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and W is a potential depending on \ an d 7. The determinant g is computed 
from 

9 = (detj^ 1 })- 1 (6.6) 



and comes out from ( |3.19|) as 



g(s; t)- 1 = - 3 {(4r 2 3 + 27r 3 2 ) + 2 S r 2 (r 2 3 + 9r 3 2 ) 

This function can be expressed by {r^} easily as follows from ( fO 

1 



9{s\t{t)) 



--(4rf + 27rf)-(detfi(r')) 2 



(6.7) 



(6.8) 



After multiplication with appropriate powers of s we obtain from ( [3~8] 



r(£,r/) = --D(£V)(dettf(£W)) 2 

where T is now the determinant of the matrix ( |4.5| ). 
In order to fix \ we start from 



e x Qe x 



^ / d dx \ _i / <9 9x N 



(6.9) 



We match the first order differential operator parts in 



(6.10) 
3|) and flBTTOl) 



- 2 7^9ab +m = -22 7T-( ln \f9)9ab 



a dr a ' - dr a 

The integrability condition for this differential equation is 

d ^ d ^ 



(6.11) 



(6.12) 



If (|6.12| ) is fulfilled we obtain a function p so that 



-1 d 
h u 'b 



In fact, we can show easily that 



P = In y/g 



(6.13) 



(6.14) 
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and 

r 2 = 3 + 2st 2 

r 3 = 2sr 3 (6.15) 
fulfill Q5TT2D and ( |5TT3| ). Then (|67TT| ) is solved by 

x = I(l + 7 )l nv ^ ( 6 .i 6 ) 
and the potential comes out as 

+TE^ (6-17) 

Due to (|6.13|) , (|6.15|) the first term in the potential is a constant and we get 

l/o x ^ i d In ,/q d In . /o 
W(r) = i ( 7 2 - g T fr. + const ( 6 - 18 ) 

a, b a " 

The wave functions 

^ = e" x ^ = g~\<W>(p (6.19) 

have a zero of order 5(1 + 7) on each trajectory C n and on the separatrix 
which can be seen from 

g-2 = cons t. \ Y[smJ-{ yi - yj )\ (6.20) 

which follows from ( |4.17| ), ( [4.24| ), ( |4.30| ) and ( |6.9| ). Moreover the spectral 
degeneracy is lifted if 7 is irrational (and s = 1, say) due to the form (|6.15|) 
of r 2 and r 3 (compare with ( |5.2(J| )). Thus we conclude that 

(1) for 7 = the Schrodinger operator is only formally selfadjoint, since 
its eigenfunctions form an incomplete set; 

(2) for 7 > the Schrodinger operator is self-adjoint if 7 ^ Q at least. For 
(see O) 

7 = 2z/-l, v G 7Z (6.21) 

insertion of ( |5.20| ) into (|6.18| ) gives the Sutherland model 0, which is 
known to be exactly soluble. 

The wavefunctions ip have support on any minimal triangle of a sector E g in 
y-space. This triangle is bounded by three trajectories 

C C C 

n\ <n 2 < n 3 , n 3 = ni + n 2 ± 1 (6.22) 
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7 The AT = 4 Calogero model 



The Calogero model for N = 3 is still too special to draw general conclusions 
from it, e.g. its Hamiltonian is separable. In fact we were able to find 
flat quadratic Lie algebraic forms also for N = 4. With s, w 3 , as in ( p,16|) , 
( |3.17| ) we make the ansatz 



922 = ~ 2r 2 - (jl + aiw 4 r 4 )s 
—(a 2 wlr^ - 03^4 T 2 r 4 ) s 2 

—a 4 w 2 r1s' i (7.1) 

923 = ^ 3r 3 - ar T 2 T 3 S - ttQWiT^S 2 (7.2) 

92.1 = - 4r 4 - (07^2X4 + agw^wlr^s 

—agw^r^s 2 (7.3) 



933 = -4w 4 w 3 2 r 4 + w 3 2 t 2 

-(a 10 w i w 3 2 T 2 T 4: + anr^s 

2 -2 2 2 i'-r a \ 

-a 12 w A w 3 r 4 s (7.4) 

9u = +]^w1 1 t 2 t z - a 13 T 3 r A s (7.5) 
3 

9ll = -2w 4 ~ 1 t 2 t 4 + -wjw^r 2 - 0,147-4 s (7.6) 

Here w 3 and W4 can be scaled to one by the automorphisms (|3.14j) , ( |3 . 1 5|) . 
Setting then s = the Calogero Riemannian for N = 4 appears. 

We have to solve the equations resulting from equating all six components 
of the curvature tensor to zero. These are polynomial equations and depend 
on the fourteen variables {a a } of (|7.1|)-(|7TBD. The number of equations lies 
somewhere between 10 3 and 10 4 . Using the seventy-eight variables {C Q( g} 
from ( |2.4| ) the number of equations would have been considerably bigger 
still. We solved the equations by solving first all equations with less than 
thirty terms and inserting then the result into the bulk of the equations. This 
trivialized almost all of them. The result is 



Oi = 


3ai — 5 




a 8 = 


-|a 7 + 


5 


a 2 = 


1„2 n , 1 

2 U 7 U7 T" 2 


11 

16 


a 9 = 


-\a 2 + 


|a 7 


o 3 = 




aio = 


a 7 — 1 




04 = 


-±a| + a 2 - 


4°7 + 2 


an = 


1 




a 5 = 


-\on + 4 




ai2 = 


-\a 2 + 


a 6 = 


gUy 2 u 7 ~r 


3 
8 


ai3 = 


|a7+| 




a-j = 


a-j 




ai4 = 


—07 + 3 





(7.7) 
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Given the Riemannian (|7.1|) - (|7.7|) and the unperturbed form 



r (o) 



-2£ -37fe 

-377; + 

. / 1 y.1 I 

-4^4 2^ % 



-4r7 4 



\ 



1 ^' ' 

-2^ + !% 2 / 



we can solve (14.81) in terms of three functions 



with the Jacobian matrix 



H 



Vs = G 3 {£', 773,774) 
7/4 = G 4 (t' ,ri 3 ,ri 4 ) 



/ 8F dGa dG 4 \ 

( I I I 

\ 9»?4 9»7 4 d?^ / 



We set 

U>3 = W4 = 1 

We find that only F depends on 07 linearly 

F = F 1 + a 7 F 2 

and that 

ro = —Ga 

2 

Consequently detif is independent of a 7 . The functions F^G^Ga 
expressed in summed up form correspondingly to ( 4.14j) ,( I.15|) 



E 

perm, of {1, . . . , 4} 



sin { 



S \ 2 



Vi] sin { 



S \ 2 



X COS { 



S \ 2 



t/aJcos{ 



S\2 



E 

perm, of {1, . . . , 4} 



sin { 



S \ 2 



Vi} s in { 



S \ 2 



2/i> 



S \ 2 



x sin { - 7/ fc } cos { 



S \ 2 



1//} 



Gi 4sin{( |V yi }sm{(~y y 2 } 



x sin { 



1/3} sin { 



2/4} 
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The separatrix 

D = (7.19) 
is given by the Vandermonde determinant 

D = -s 6 V( yi ,y 2 ,y 3 ,y,) 2 (7.20) 

or by 

D = -4detr (0) (7.21) 
D = - 256r]l + 128£ 2 r^ 

-16CV + 4f% 2 - 144^ 3 V (7.22) 



Moreover we obtain 



with 



detH = [] (7-23) 



l<j<j<4 a *J 



«0-= ( TiYiVi-Vi) ( 7 -24) 



and 



detr(e, 773,^4) = -\D(Z\Th,V4)<teH(Z',Th,ritf (7.25) 

In the case N = 4 the discussion of the Schrodinger operator proceeds along 
the lines of section 5. The linear differential operator part in ( |6.1|) is now 
(for s = w 3 = W4 = 1) 

r 2 = i(a 7 + 5)r 2 -i(a2-4a 7 + 3)r 4 + 6 (7.26) 
r 3 = 3r 3 (7.27) 
r 4 = r 2 -i(a 7 -9)r 4 (7.28) 

The integrability condition (|6.12|) is fulfilled and the solution of the differ- 
ential equation ( |6.13| ) is ( |6.14j ) as for N = 3. The potential formula ( |6.18| ) 
remains correct and yields the N = 4 Sutherland model potential. 



8 Concluding remarks 

In carrying over these results to general iV we must clarify the meaning of 
the parameter a 7 in (|7.1|)— (|7T7[) . Assume 

Vk = G k (£',r)' 3 ,...,r]' N ) (8.1) 
k e {3,4,..., N} 



22 



to be a basic diffeomorphism satisfying the equation analogous to (|4.8|) . Then 
another basic diffeomorphism is obtained from (O) by application of a linear 
map with triangular matrix 



N 

£ = F + ^ 0,2 m G m 
m=4 
m even 
N 

Vk = Gk+ O-kmGm (8-2) 

m=k+2 
m-k even 

So the flat Riemannian obtained by our algorithm depends on -r{N — 2) 2 
(N even) respectively 4(iV — 1)(N — 3) (AT odd) parameters, which drop out 
of det H and thereby drop out of the Schrodinger operator. Despite of the 
trivial appearance of these parameters in the diffeomorphism fl3.2|) , they are 
connected with an annoying complication of the Riemannian as (7.7) shows. 

We have not discovered yet new exactly solvable models nor have we 
been able to transform new known integrable models into exactly solvable 
forms. However, the exactly solvable form of the Sutherland models found 
here differs from the known one ||. Moreover, we have shown that the 
formalism of the flat quadratic Lie algebraic forms implies an algorithm for 
the deformation of exactly solvable models into other exactly solvable models. 
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